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The Gene Ontology (GO) provides biologists with a controlled terminology that describes how 
genes are associated with functions and how functional terms are related to each other. These 
term-term relationships encode how scientists conceive the organization of biological functions, and 
they take the form of a directed acyclic graph (DAG). Here, we propose that the network structure 
of gene-term annotations made using GO can be employed to establish an alternate natural way to 
group the functional terms which is different from the hierarchical structure established in the GO 
DAG. Instead of relying on an externally defined organization for biological functions, our method 
connects biological functions together if they are performed by the same genes, as indicated in a 
compendium of gene annotation data from numerous different biological experiments. Grouping 
terms by this alternate scheme provides a new framework with which to describe and predict the 
functions of experimentally identified sets of genes. 



1. INTRODUCTION 

The Gene Ontology (GO) 0][29] has been around for 
over a decade, during which time it has been widely uti- 
lized both to validate and to predict the results of biolog- 
ical experiments (see, for example [9j El EE [HI EJ [33j 
[34]). The structure of the ontology, where different "cate- 
gories" or terms are related to each other in a hierarchical 
fashion, provides a well-established format with which to 
classify and subclassify all biological functions and pro- 
cesses. This classification approach is well-structured and 
well-characterized, however, we seek to determine if it is 
the only natural way in which to classify this type of 
biological information, or if any other, alternate natural 
groupings of functions might exist. We address two main 
questions. First, does there exist another natural way to 
organize the functional terms that is distinct from the 
ontological organization? Secondly, if such an alternate 
classification exists, can it be used to interpret biological 
data? 

Our approach is outlined in Figure 1. We begin by 
considering the term relationships as defined by the GO 
hierarchy. We then add in annotation information that 
reflects gene-term relationships obtained from numerous 
different biological experiments. We encapsulate these 
connections in the form of a bipartite network. Next, 
we use this bipartite network to construct another net- 
work describing the relationships between the functional 
terms based on shared gene annotations. We then apply 
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standard community structure finding algorithms to par- 
tition this annotation-driven network into communities of 
terms. We compare these communities to the previously 
defined branches from the GO hierarchy and show that, 
although there are some similarities, there are also very 
strong differences between the two ways of organizing 
functional terms. Finally, using functional enrichment 
techniques, we find a strong association between both 
our functional term communities and GO branches (on- 
tological groupings of terms) with experimentally-derived 
sets of genes. 



In the past there has only been minimal investiga- 
tion of how biological functions might be related to each 
other outside of the established ontology structure, and 
the majority of this has focused on discovering individ- 
ual links between functions [16] rather than investigating 
the structure as a whole. In this work, we propose an 
alternate, viable classification of functional terms that 
relies on the network structure of gene-term annotations 
rather than ontological relationships. We test the appli- 
cability of this classification, utilizing functional analysis 
techniques in order to evaluate the enrichment of cancer 
signatures (sets of genes associated with a particular can- 
cer) in both term communities and GO branches. We find 
that certain signatures are highly enriched in the found 
communities and not GO branches. We therefore sug- 
gest that by linking GO functional terms based on shared 
genes, we can create an alternate, biologically meaning- 
ful, network- derived organization of the functional terms 
that is both distinct from the established GO DAG and 
can also be used to investigate biological systems. 
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FIG. 1: Visual representation of our approach. First, we sum- 
marize gene annotations made to functional terms in the Gene 
Ontology Hierarchy as a gene-term bipartite graph. From 
these gene- term relationships, we create a projected term- 
term network. We partition the term-term network into com- 
munities and compare those term communities to branches of 
terms in the DAG. Finally, we perform functional enrichment 
analysis on gene sets using both these term communities and 
GO branches. 



2.1. 



2. BACKGROUND 



Identifying and Comparing Community 
Structure in Networks 



In recent years, complex networks tools have been 
used alongside traditional bioinformatics techniques to 
study many different kinds of biological networks [24] . 
including, but not limited to, gene regulatory networks 
[2QJ[28], protein-protein interaction networks p~5j[3l], and 
metabolic networks [13] [35]. Developments in network 
theory provide the computational tools needed to calcu- 
late global properties of such networks, lending insights 
into the behavior of the systems represented by these 
networks. Many networks exhibit community structure, 
meaning that there are clusters of nodes in the network 
within which there are many edges but between which 
there are few edges. Recent developments in network 
theory have been able to detect such functional modules 
[23] in a computationally efficient and accurate manner 
[6]. Algorithms for detecting network structure, although 
most fully developed for un-weighted, un-directed graphs, 
can also be applied to weighted [22] or directed [18] net- 
works with only slight modifications. 

In order to quantify the strength of community struc- 
ture we use a quantity known as modularity [22] [23] . 
Modularity (Q) can be denned as: 
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where S is the Kronecker delta function, x* is the com- 



munity of node z, ki is the degree of node z, A is the 
adjacency matrix, a matrix with values representing the 
weight between nodes i and j, and m is the total weight 
of the edges in the network [3 . Traditionally, in order to 
partition the network into communities, the resolution 
parameter, r in Equation [T] is set equal to zero. Vary- 
ing this value allows one to look for alternate divisions 
of a network into communities at different scales, with 
r > uncovering sub-structures in the network, finding 
new communities, or breaking apart communities found 
at lower resolutions. 



2.2. The Gene Ontology 

Ontologies are utilized across many disciplines includ- 
ing economics, artificial intelligence, engineering, library 
science, and biomedical informatics (for a few examples, 
see [5j [19l [27]). They are a structural framework for 
organizing information, representing knowledge and de- 
scribing the relationship between different ideas. Fur- 
ther, they provide a method for conceptualizing common 
features shared across multiple cases. The Gene Ontol- 
ogy, specifically, describes the relationships between dif- 
ferent biological concepts or functions [4] . It breaks these 
concepts into three main domains, or distinct ontologies: 
"Biological Process" (BP), describing sets of molecular 
events, "Molecular Function" (MF), describing the activ- 
ities of gene products, and "Cellular Component" (CC), 
describing parts of a cell or its external environment. 
Each of the three primary domains in GO takes the form 
of a directed acyclic graph (DAG), in which "child" func- 
tional categories, or "terms" , are subclassified under one 
or more "parent" terms, using "is a" and "part of" re- 
lationships. Each parent and all its subsequent progeny 
therefore define multiple, overlapping, sets of terms, or 
"branches" in GO. Note that a child term classified in one 
of the primary domains cannot have a parent classified 
in a different primary domain. Further the primary do- 
mains vary greatly in size with the BP domain containing 
approximately two-thirds of all terms. Using GO, genes 
are annotated to terms representing their particular role 
in a cell, and these annotations are transitive up the re- 
lationships in the DAG such that each "parent" term 
takes on all the gene annotations associated with any of 
its progeny [30] . 

In the following analysis we explore if there exists an- 
other, natural way to classify terms independent from the 
ontology structure. To begin, we use term-term ontol- 
ogy relationships and gene-term annotation information 
for human genes downloaded from the GO website (ge- 
neontology.org). We determine the terms with at least 
one human gene annotation, and plot the cumulative 
distribution of the size of the branches corresponding 
to each of these terms and its annotated-progeny (Fig- 
ure 2(a)). The heavy-tailed distribution is a result of 
the hierarchical DAG structure as the members of each 
branch are also members of the parent branch (es). Be- 
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FIG. 2: A comparison of branches in the GO DAG and term 
communities found by partitioning the term network. (A) 
The cumulative Distribution for the sizes of all branches in 
the Gene Ontology and all unique term communities found at 
the various resolutions. These distributions demonstrate that 
the number and sizes of branches and communities is similar. 
(B) Distribution of J m , the maximum similarity a commu- 
nity or branch with ten or more members has compared to 
all other branches or communities with ten or more members, 
respectively. Although a small number of communities and 
branches have similar memberships, most are highly dissimi- 
lar. 



cause of the transitive nature of gene annotations, larger 
branches, which generally represent broader terms, often 
have many genes annotated to them, whereas terms with 
fewer or no progeny generally represent more specific bi- 
ological activities that are performed by relatively fewer 
genes. Therefore, as with the branch size, a distribution 
of the number of genes annotated to each term follows a 
heavy-tailed relationship; however, the number of terms 
to which each gene is annotated does not appear heavy- 
tailed (see [Hj [12]). 

Next, we used the annotation information to construct 
a gene-term bipartite network. We represent this network 
in the form of an no x Ut adjacency matrix, where tlq is 
the total number of genes and tit is the total number of 
terms listed in the annotation file. In this matrix a value 
of one indicates a known connection between the corre- 
sponding gene and term, and a value of zero indicates 
that the gene is not associated with that term. Thus, 

J 1 if gene p is annotated to term i 

JDpi — \ . (2) 

I if gene p is not annotated to term i 

This bipartite network represents a human-specific sum- 
mary of the relationships between genes and terms. 

We note that although GO is broken into three primary 
domains and gene-annotations are made to the ontology 
for many species, for simplicity in the following analysis 
we will combine information from all three domains and 
use annotation information only that pertains to human 
genes. 



3. METHODS 

3.1. Projecting Term Networks based on Gene 
Ontology annotations 

In this section we use gene-term annotations to con- 
struct a network representing term-term relationships. 
Using the bipartite network (Equation |2| one might cre- 
ate a term network by simply joining together any pair 
of terms that share common genes; however, this ap- 
proach would lose a large amount of information as con- 
nections between high degree terms would be given the 
same weight as connections between low degree terms. 
We correct for the skewed term degree distribution by 
constructing a diagonal weighting matrix, with elements: 

= (3) 

q=l 

or simply one over the degree of the term in the bipartite 
graph. Using w, we can project a term- network, T, whose 
edges are modified by this weighting matrix: 

no 
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In this network the weights of T have a maximum value of 
one when the same single gene is annotated to both term 
i and term j and a minimum value of zero when none 
of the genes annotated to term i are annotated to term 
j. The use of the weighting matrix biases the weights 
of network edges to those between low degree terms and 
therefore the lower branches of the DAG. 



3.2. Identifying Communities of GO terms 

Our first goal is to identify the community structure in 
annotation- driven term-term relationships, meaning we 
wish to find clusters of terms within which there are many 
or high-weight relationships in our projected network, 
but between which there are only few or low- weight rela- 
tionships. To determine the community structure of the 
network we used a weighted version of the Fast Greedy 
Community Structure algorithm [6] and determined the 
communities of terms at maximum modularity. Fifty- 
six communities were found with a modularity value of 
Q = 0.57, indicating that there are indeed many individ- 
ual modules of terms which share gene annotations. 

The branches in GO represent multiple, overlapping 
sets of functional categories rather than one discreet par- 
tition of terms, indicative of functional structure at many 
different levels of specificity. Since standard community 
structure approaches only find one partition of network 
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FIG. 3: Visualization of communities (circles) of GO terms found at the six lowest levels of resolution (rows), in increasing 
order (top to bottom). Width of line connecting two communities is proportional to the percentage of terms in the child 
community that are also in the parent community. The size of communities is proportional to the log of the number of terms 
in the community. Color represents the normalized percentage of terms in the community which belong to the BP (yellow), 
MF (magenta) and CC (cyan) primary domains. 



nodes (in our case terms), we implemented a modified 
version of the Fast Greedy that maximizes modularity 
for non-zero values of the resolution parameter (see Equa- 
tion [I]) in order to find many different viable partitions. 
We varied this parameter several orders of magnitude to 
find 11491 different communities. The different values 
of the resolution parameter were chosen to give commu- 
nity sizes that were roughly similar to those defined by 
the branches at different levels of the GO DAG. Like GO 
branches at different levels, term communities at different 
resolutions are highly overlapping. We gave our commu- 
nities numeric identities that vary from TC: 0000001 to 
TC:0011491 and will refer to them as such in the follow- 
ing analysis. The cumulative distribution of the number 
of members in these communities is shown in Figure 2(a). 
As with the GO branches, this is a heavy-tailed distribu- 
tion. 



RESULTS: AN ALTERNATE "NATURAL" 
GROUPING OF GO TERMS 



communities. 

While partitions at low resolutions are somewhat simi- 
lar, we find that high resolution partitions are vastly dif- 
ferent. Communities at higher resolutions do not merely 
represent the "splitting apart" of communities at lower 
resolutions (represented by a child community only con- 
necting to a single parent), but instead each resolution 
often brings about a new way of partitioning the network. 
Sometimes members from multiple communities found at 
lower resolutions split and/or combine with other lower- 
resolution communities to form a new community at a 
higher resolution. 

The colors of the communities illustrate how the struc- 
ture of annotation-driven term relationships is distinct 
from the structure of those relationships as defined by 
GO branches. Each GO branch can only belong to one 
primary ontology, and thus would be pure cyan, magenta 
or yellow in this type of visualization, however, commu- 
nities, even smaller ones and those found at higher reso- 
lutions, generally contain members from multiple ontolo- 
gies, resulting in a rainbow of colors. 



4.1. Illustrating Community Structure at different 
resolutions 



To better understand the relationships between the 
communities found at different resolutions, we visualized 
the term communities with ten or more members for the 
six lowest values of resolution used (Figure 3). In this 
visualization each community is represented by a single 
circle, whose radius scales as the log of the number of 
terms belonging to that community and whose color cor- 
responds to the percentage of members from each pri- 
mary domain that belong to that community. Between 
the communities found at adjacent resolutions, we draw a 
line from a community at a higher resolution to a commu- 
nity at a lower resolution if at least 10% of the members 
of the community from the higher resolution also belong 
to the community at the lower resolution. The thickness 
of the line is indicative of the overlap between the two 



4.2. Comparing Term Communities and GO 
Branches 



Next we compared the term communities with 
branches of the GO DAG. In order to quantify the spe- 
cific differences between the communities and branches, 
we directly compared the membership of each commu- 
nity with the membership of all branches in GO using 
the Jaccard similarity (J), calculating a similarity value 
for every community-branch pair that takes the value: 



J(x,y) 
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For each community (x), we determined the correspond- 
ing branch (y) that has the highest overlap in member- 
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ship by this measure: 



J m (x) = max J(x,y), 



(6) 



and vice versus. 



Because the exact value of the Jaccard similarity is 
highly sensitive to incremental changes in set member- 
ship when comparing sets with only a few members, we 
will limit all the following analysis to communities and 
branches that contain ten or more terms in order to focus 
on the most robust results. Figure 2(b) shows the dis- 
tribution of J m comparing these 2370 communities and 
2151 branches. Although a handful of communities and 
branches are quite similar to each other, the majority 
of communities are dissimilar to the GO Branches and 
vice versus. To better interpret these values, we selected 
several communities to inspect more closely. 

First we selected a community with a very high J m 
value to inspect (Figure 4(a)). TC:0007391 is most sim- 
ilar to GO:0070570 ("regulation of neuron projection re- 
generation") with J m = 0.6667. It is interesting that in 
addition to members from the BP domain TC:0007391 
also includes two members from the MF and CC do- 
mains, "neutrophin receptor activity" and "perineuronal 
net" respectively, the former of which is involved in the 
regeneration of injured axons [10 while the degradation 
of the latter has been shown to favor axon regeneration 
[26] . This indicates that terms found in the community 
but not the branch are consistent with known biology. 

Next we selected TC:0001295 and GO:0052547 to il- 
lustrate the shared information typically found between 
a community and the branch (Figure 4(b)). The branch 
defined by GO:0052547 has members that belong to seven 
distinct communities, demonstrating that not only are 
communities often distinct from branches, within the 
branches themselves the annotation-driven classification 
is often very distinct from the defined ontological rela- 
tionships. One might suppose that a significant factor 
to the dissimilarity found between TC:0001295 and the 
branches in GO is partially attributable to the fact that 
it has term members that belong to all of the primary do- 
mains. To address this issue for our last example we illus- 
trate a term community (TC: 0000936) whose members 
all belong to the same primary domain (Figure 4(c)). As 
with TC:0001295 this community only shares a handful 
of members with its corresponding branch (GO:0060538). 

From this analysis we conclude that although there is 
occasional similarity between our found communities and 
GO branches, the communities are not simply a recapit- 
ulation of the DAG. We have repeated portions of this 
analysis constructing the term network and correspond- 
ing partitions three more times, using only annotations 
specific to only one of the three primary ontologies, and 
observe similar results (data not shown). 
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(a)TC:0007391, J m = 0.6667 with GO:0070570 





(b)TC:0001295, J m = 0.1818 with GO:0052547 





(c)TC:0000936, J m = 0.1667 with GO:0060538 



FIG. 4: Three example comparisons between communities 
and branches: (A) TC:0007391 compared to GO:0070570, 
(B) TC:0001295 compared to GO:0052547, (C) TC:0000936 
compared to GO:0060538. In each panel on the left hand 
side a community and its inter-community connections in the 
annotation-driven term network is shown and on the right 
hand side the branch with which that community has the the 
highest Jaccard similarity is illustrated. In the right panel 
edges represent the ontological associations defined by the 
Gene Ontology term hierarchy. Each term member of the 
community or branch is colored both by its associated primary 
domain (inner color - BP:yellow, MF:magenta, CCxyan) and 
its community membership (outer color), determined at the 
same resolution value as the illustrated community. Terms 
common between the community and the branch are circled. 



4.3. Capturing the Biological Information in Term 
Communities 



We have illustrated that our term communities rep- 
resent a natural partitioning of biological functions that 
is distinct from the GO DAG, however, the biological 
meaning of these communities is, at this point, unclear. 
On a mathematical level they represent sets of biologi- 
cal functions that are generally performed by the same 
collection of genes. Labeling and understanding the bio- 
logical meaning behind these communities is vital if they 
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are to have the same of wide-range applications as the 
GO branches. 

At first one might consider determining the infor- 
mation contained within each community by manually 
inspecting the functional categories belonging to those 
communities. For communities containing only a hand- 
ful of terms, such as those illustrated in Figure 4, this is 
doable. Unfortunately this task can be daunting at best 
for larger communities, each of which contains several 
hundred members. Therefore, in order to easily interpret 
the contents of our communities we summarize the in- 
formation contained in each in the form of word clouds 
(Figure 5(a)). Specifically, for each community, we deter- 
mine the descriptions corresponding to the member terms 
of that community. We then count the number of times 
an individual word appears across all these descriptions 
and calculate the statistical enrichment of the frequency 
of that word in the community compared to its frequency 
across the descriptions of all GO terms. Finally, we use a 
publicly available word-cloud making program [Tj to illus- 
trate the words belonging to each community, scaling the 
relative sizes of the words based on the word's statisti- 
cal enrichment in the community and coloring each word 
based on the normalized percentage of times the term is 
it derived from is a member of each primary domain. 

As an example, take Community TC:0000400. Looking 
at all 335 members of this community and trying to man- 
ually determine what biological information they repre- 
sent could be difficult or impossible. The word cloud pre- 
sentation (Figure 5(b)), however, easily summarizes this 
information and reveals that this community includes bi- 
ological concepts related to various types of RNA, includ- 
ing "rRN A" , "tRNA" , "mRN A" , "LSU-rRNA" , "SSU- 
rRN A" , "ncRNA", "RNA-polymerase" and more. It has 
members belonging to all three primary domains, and 
often the individual words shared between those mem- 
bers are associated with multiple domains, resulting in 
a complex coloration. We note that this community is 
not highly similar to any particular branch in GO. It has 
the highest similarity to GO:0016070 or "RNA metabolic 
process", with a value of J m — 0.22146. Other word 
clouds show a similar richness of information. For ex- 
ample TC: 0000061 contains many words related to the 
heart such as "cardiac", "muscle", "ventricle", "ventric- 
ular" and "heart." (Figure 5(c)). This community is also 
distinct from the branches with a J m value of 0.149 with 
GO:0072358. 

One can also represent the biological information con- 
tained in branches in the form of word clouds, although, 
because the members of each branch can only belong to 
one of the three primary domains, all the words in the 
cloud will be the same color. Two branches (GO:0000003 
and GO:0002376) are illustrated in Figure 5(d)-(e). The 
first clearly contains terms pertaining to sex-related pro- 
cesses as it contains words such as "female", "sex", 
"prostate" and "male." This corresponds well to the 
name of the parent term in the branch, "reproduction". 
Similarly, the cloud for GO:0002376, whose parent term 
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(d)GO:0000003 Word 
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FIG. 5: (A) A schematic illustrating how the biological infor- 
mation contained in a particular community or branch can 
be summarized in a word cloud. (B-E) Term Communi- 
ties (TC:0000400, TC:0000061) and branches (GO:0009607, 
GO: 0050896) summarized as word clouds. In each case the 
color of a word represents how often the term description 
containing that word belongs to each of the primary d oma ins 
(BP: yellow, MF: magenta, CCxyan, also see Figure 2.2 for 
mixed-domain coloration) and size represents that word's sta- 
tistical enrichment in that community /branch. 



name is "immune system process" contains words per- 
taining to the immune system. 



4.4. Using the term communities to evaluate and 
predict genetic function 

Finally, we wanted to test how our communities might 
be used in one common application of the Gene Ontol- 
ogy: functional enrichment analysis. To perform this 
analysis we chose to use the structure-preserving form of 
Annotation Enrichment Analysis (SP-AEA) [IT] since it 
has been shown to better estimate the biological func- 
tions of experimentally derived sets of genes, and has a 
conceptual framework conducive to estimating functional 
enrichment between sets of terms and genes, rather than 
simply between two sets of genes. 

To begin, we downloaded a collection of experimentally 
derived genes sets from the Gene Signatures Database 
(GeneSigDB) [8 j. This database is a manual curation of 
previously published gene expression signatures, focusing 



7 



Ovarian (Ouellet '05) 
Bone (Heller '08) 
Sarcoma (L( 
Leukemia (Walker '06) 
StemCell (Osada '05) 
HeadandNeck (Bredel '06) 
Sarcoma (Missiagli 

Viral (Urosevic '07) 
EmbryonicStemCell (Xu '09) 
Bladder (Os 
Sarcoma (Filion '09) 
Breast (Creight( '"" 
Intestine (Vecchi '07) 
Breast (Mazurek '05) 
HeadandNeck (Roepman '06) 
Lung (Hou '10) 
Breast (Mira '09) 
Colon (Grade '06) 
Leukemia (Calln 



1 ■ 


m 


■ i 








■ 1 


r 


■ 


■ 


■ 

■ ■ 


r 










■ 1 



n 



Sm?^™o8ooooooo™Sooooooooo8oo 

ooo-*co-*oldoooooooooooooooooooooooooooooooooo 

8888SSSS88888888SSSSS88888S88888888888 



FIG. 6: A heat map showing the statistical enrichment of 
selected cancer signatures (see text) in GO branches and term 
communities. 



anient" (TC:0000249), "morphogenesis" and "develop- 
ment" (TC:0000365), "blood", "pressure" and "contrac- 
tion" (TC:0000582), with the other clouds generally con- 
taining these words in different combinations (see for ex- 
ample TC:0000061, illustrated in Figure 5(c)). 

The second signature is a list of bladder cancer spe- 
cific genes [25]. Most of the words emphasized by the 
community clouds are related to cell proliferation. For 
example, it is enriched in TC:0000400 (illustrated in Fig- 
ure 5(b)) and the other clouds emphasize words such 
as "cell-cycle", "mitotic", "meiotic", "checkpoint", "re- 
pair", "replication", "recombination", "telomere", "spin- 
dle" , "complex" , "DNA" , "chromosome" , "histone" , and 
"methylation" . Although the connection to the bladder 
is not obvious, the connection to cancer and the high rate 
of cell proliferation in tumor cells [2] is apparent. 



primarily on cancer and stem cell signatures [7]. In the 
following analysis we will use all 509 human signatures 
from this database that contain at least 100 and less than 
1000 genes annotated in the Gene Ontology. Using SP- 
AEA (with one million randomizations), we estimated 
the significance of enrichment for these 509 signatures in 
both the term communities and GO branches. 

Term communities show a level of statistical enrich- 
ment equal or greater to that of GO branches, lending 
biological validity to our term communities. Furthermore 
this level of enrichment is not evident using random com- 
munities (data not shown). We wished to know if there 
was a context in which our term communities captured 
biological information missed by the branches, or vice 
versus. Thus we selected gene signatures that were sig- 
nificantly enriched (p < 10 -6 ) in at least one commu- 
nity/branch but not significantly enriched (p > 5 x 10 -5 ) 
in any GO branch/community. Figure 6 shows a heat 
map of the nineteen signatures that met this criteria 
with the communities and branches for which they are 
enriched. 

It is immediately striking that of these signatures, 
the majority are enriched in communities and not GO 
branches, rather than the other way around. Of course, 
some of these communities may capture redundant in- 
formation and only vary by a few members, however, 
we point out this is equivalent in structure to the GO 
branches, where the branch of a term will often contain 
virtually the same membership as that of its immediate 
progeny and parent (see, for example, the branch defined 
by GO:0060538, illustrated in Figure 4(c)). 

Two signatures, in particular, are enriched in a col- 
lection of communities. The first, an embryonic stem 
cell signature [32 , represents genes that are up-regulated 
in cardiomyocytes compared to non-selected embryoid 
bodies and hESC. The communities represented in this 
signature contain several different themes, all consis- 
tent with the expected properties of genes selected 
from stem cells and related to the heart. The cor- 
responding clouds emphasize words such as "cardiac" 
and "muscle" (TC:0000012), "actin", "myosin", and "fil- 



5. DISCUSSION 

The network structure of gene annotations using the 
Gene Ontology has not previously been exploited in a 
manner that reveals an organization of biological func- 
tion that is unique from the published hierarchical clas- 
sification of the GO DAG. By analyzing functional anno- 
tation data we were able to construct an alternate, natu- 
ral, and biologically-relevant way in which to categorize 
cellular functions. This categorization is structurally and 
conceptually distinct from the GO DAG and allows for 
functional relationships between terms that do not share 
a parent/child relationship. It takes advantage of a large 
amount of data from a variety of sources and creates a 
classification scheme that is motivated primarily by the 
data reported rather than the organization of human con- 
ceptions. 

The term communities denned in this work represent 
an integration of information across all three primary do- 
mains in GO that, to the authors' knowledge, has not 
previously been investigated systematically. Using the 
simple principle of co-annotation we suggest that in the 
future biological concepts from other smaller databases 
could also be analyzed or even combined with these re- 
sults. Additionally, we concede that the communities de- 
fined here likely do not represent the only way to group 
functional terms outside of the ontology structure. Even 
so, we believe that the functional enrichment analysis 
presented here demonstrates that the term communities 
we define are clearly more than a mathematical phe- 
nomenon and have a high potential to be used to better 
interpret biological data. 
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